Domain State Model for Exchange Bias 



U. Nowak, A. Misra, and K. D. Usadel 
Theoretische Tieftemperaturphysik, Gerhard-Mercator-Universitat Duisburg, 47048 Duisburg, 

Germany 
(February 1, 2008) 

Abstract 



Monte Carlo simulations of a system consisting of a ferromagnetic layer 
exchange coupled to a diluted antiferromagnetic layer described by a classical 
spin model show a strong dependence of the exchange bias on the degree 
of dilution in agreement with recent experimental observations on Co/CoO 
bilayers. These simulations reveal that diluting the antiferromagnet leads 
to the formation of domains in the volume of the antiferromagnet carrying 
a remanent surplus magnetization which causes and controls exchange bias. 
To further support this domain state model for exchange bias we study in 
the present paper the dependence of the bias field on the thickness of the 
antiferromagnetic layer. It is shown that the bias field strongly increases 
with increasing film thickness and eventually goes over a maximum before 
it levels out for large thicknesses. These findings are in full agreement with 
experiments. 
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For a ferromagnet (FM) in contact with an antiferromagnet (AFM) a shift of the hystere- 
sis loop along the magnetic field axis can occur which is called exchange bias (EB). Usually, 
EB is observed after cooling the system with the FM magnetized in saturation below the 
Neel temperature Tn of the AFM or after cooling the entire system in an external magnetic 
field. Although this effect is well known since many years |I],|2| its microscopic origin is still 
discussed controversially. For a review see a recent article by Nogues and Schuller ||. 

In a previous Letter Q we have reported on EB observed experimentally in Co/CoO 
bilayers as a function of volume defects in the antiferromagnet. Of particular importance in 
this study was the observation that it is possible to strongly influence EB in Co/CoO bilayers 
by diluting the antiferromagnetic CoO layer, i. e. by inserting nonmagnetic substitutions 
(Coi-^Mga-O) or defects (Coi-^O) not at the FM/AFM interface, but rather throughout the 
volume part of the AFM. While the undiluted samples show only a very small EB, dilution 
increases EB dramatically. Since for all samples investigated a 0.4 nm thick CoO layer with 
minimum defect concentration was placed at the interface the observed EB is primarily not 
due to disorder or defects at the interface. Rather, the full antiferromagnetic layer must 
be involved and we have argued that in our systems EB has its origin in a domain state in 
the volume part of the AFM which triggers the spin arrangement and thus the FM/AFM 
exchange interaction at the interface. Indeed, in diluted antiferromagnets when cooled in 
external fields metastable domains occur carrying a surplus magnetization and having a 
very slow dynamics (for reviews see |5|||). These domains are frozen to a large extend 
during hysteresis cycles and their frozen magnetization is the origin of EB. This domain 
state model for EB was supported by large scale Monte Carlo simulations performed at 
finite temperatures. ||,[7| 

To gain additional evidence for the domain state model in the present paper we will 
concentrate on the dependence of EB on the thickness of the AFM. The system consists 
of a FM monolayer exchange coupled to a diluted AFM layer with t monolayers. The FM 
is described by a classical Heisenberg model with vector spins S_ { and exchange constant 
Jfm- The dipolar interaction is approximated by an additional anisotropy term (anisotropy 
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constant d x = — O.IJfm) which mimics the shape anisotropy leading to a magnetization 
which is preferentially in the y — z-plane. Also, we introduce an easy axis in the FM (z- 
axis, anisotropy constant d z = O.IJfm) in order to obtain well defined hysteresis loops. d z 
sets the Stoner-Wohlfarth limit of the coercive field, i. e. the zero temperature limit for 
magnetization reversal by coherent rotation (i?sw — 2d z , in our units, for a field parallel 
to the easy axis). In view of the rather strong anisotropy in CoO we assume an Ising 
Hamiltonian for the AFM. Thus the Hamiltonian of our system is given by 

H = - Jfm ±U ' £j - Yl (dzS 2 iz + d x Si X + S_i- 

— JaFM X! e i e j (T i' J j ~ ^ BejCTj 

{i,3') i 

— Jint X! e i a iSjz- (1) 

with the antiferromagnetic nearest-neighbor exchange constant Jafm < and the effective 
in-plane magnetic field B_ = Bz_+B y y. The values of the magnetic moments are incorporated 
in B and d, respectively, so that the quantities ^ denote unit vectors and = ±1 Ising 
spin variables. A fraction p of the sites of the lattice is left without a spin (quenched 
disorder: q = 0,1). For the exchange constant of the AFM which mainly determines its 
Neel temperature (also depending on the dilution, of course) we set Jafm = — Jfm/2. There 
seems to be some evidence that the exchange coupling between Co and CoO is ferromagnetic 
0] but its strength is not known experimentally. Therefore, we assume in our simulations a 
ferromagnetic coupling with (Jint = ~ Jafm)- 

We use Monte Carlo methods with a heat-bath algorithm and single-spin flip methods 
for the simulation of the model explained above. The trial step of the spin update is a small 
variation around the initial spin for the Heisenberg spins and — as usual - a spin flip for the 
Ising spins 0. We perform typically 40000 Monte Carlo steps per spin (MCS) for a complete 
hysteresis loop. To observe the domain structure of the AFM we have to guarantee that 
typical length scales of the domain structure fit into our system. For the parameter values 
used in this simulation this is achieved for systems of lateral extension 128 x 128. 
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In the simulations the system is cooled from above to below the ordering temperature 
of the AFM in an applied external cooling field B_ = B c z with B c = 0.25Jfm • The FM 
is then long-range ordered and its magnetization is practically homogeneous resulting in a 
nearly constant exchange field for the AFM monolayer at the interface. When the desired 
final temperature is reached a magnetic field B_ = Bz_ + B y y is applied which also has a 
small, constant perpendicular field component B y = 0.05Jfm in order to define a certain 
path for the rotation of the magnetization during field reversal and to avoid the system to 
be trapped in a metastable state. The z component of the field B_ is then reduced in steps 
of AB = 0.004Jfm from B = 0.2Jfm down to —B and afterwards raised again up to the 
initial value. 

Typical hysteresis loops are depicted in Fig. |l|. Shown are results for the magnetization 
of the FM (upper figure) as well as that of the AFM interface monolayer (lower figure) 
for different thicknesses t of the AFM. The hysteresis loops of the FM clearly show EB 
depending on the thickness of the AFM. The magnetization curves of the interface layer of 
the AFM are shifted upwards due to the fact that after field cooling the AFM is in a domain 
state with a surplus magnetization. This layer experienced during cooling the external field 
in addition to the exchange field of the FM both having the same direction. This shifted 
interface magnetization of the AFM acts as an additional effective field on the FM, resulting 
in EB. The hysteresis curve of the AFM interface layer follows that of the FM layer with 
a much lower saturation magnetization, however. With increasing thickness of the AFM 
layer the area of the hysteresis loop of the interface layer which is proportional to the energy 
losses in the AFM decreases indicating that the spin structure in the AFM is stabilized. 

The domain structure in the AFM interface layer is shown in Fig. |^ for an AFM consisting 
of one monolayer (upper figure) and for ten monolayers (lower figure), respectively. The 
fractal structure of these domains is obvious. It has been observed previously in bulk systems 
and was analyzed in detail |10|,|TTJ. The domain structures shown in Fig. |2| are to a large 
extend frozen. But during field cycles small spin arrangements in the domain boundaries 
can take place even at low temperatures resulting in the hysteresis loops shown in Fig. |l|. 
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The structure of the domains depends on the thickness of the AFM. For an AFM mono- 
layer the effective field acting on all AFM spins is the superposition of the strong exchange 
field and the external field. But it is well known that the size of the domains depends 
on the strength of the effective magnetic field. Large fields imply small domains and vice 
versa |T2|-|n],|llIO] . The small domains seen in the upper part of Fig. □ are thus due to the 



strong fields acting on the AFM monolayer. On the other hand if the AFM consists of ten 
layers nine of them are only exposed to the weak external field with the tendency to form 
larger domain sizes. The coupling of these layers to the AFM interface layer then results 
in a coarsening of the domains at the interface as seen in the lower part of Fig. [| Note 
that the distribution of vacancies in the interface layer is exactly the same in both parts of 
Fig. |2|. A further obvious consequence of this explanation is that the domain size becomes 
layer dependent and increases with increasing distance from the AFM interface. But after 
a certain distance from the interface the domain structure should become independent of 
the interface layer which means that also the bias field should become independent of the 
thickness of the AFM for large t. This behavior indeed is observed in our simulations. In 
Fig. ^] we show the dependence of the bias field on t for different dilutions of the AFM. The 
bias field is determined as -Beb — {B + + B~~)/2 where B + and B~ are those fields of the 
hysteresis loop branches for increasing and decreasing field where the easy axis component 
of the magnetization of the FM becomes zero. The absolute value of the bias field increases 
rapidly with film thickness t, goes eventually over a maximum and then levels out. This is 
in agreement with experiments ||. Note that for the system with the smallest dilution the 
absolute value of the bias field decreases for t > 1 much stronger with increasing thickness 
than for the more diluted films. The reason is that the less diluted systems have a stronger 
tendency to order antiferromagnetically thus reducing the net magnetization at the AFM 
interface. Again, agreement with experiments is obtained. 

In conclusion, we have found by Monte Carlo simulations further support for our domain 
state model for exchange bias. So far all our numerical results have been obtained for systems 
where the AFM has a very strong anisotropy, i.e. behaves Ising-like. It is of great interest to 



relax this condition, i.e. to consider vector spin models for the AFM. Work in this respect 
is in progress. 
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FIGURES 

FIG. 1. Simulated hysteresis loops of the model explained in the text for p = 0.4 and 
ksT = O.IJfm- The field during cooling was 0.25Jfm- Shown is the net magnetization in units of 
the saturation magnetization of the FM (upper figure) and of the interface monolayer of the AFM 
(lower figure) for different thicknesses t of the AFM. 

FIG. 2. Frozen domain states of the AFM. Shown are staggered spin configurations (grey and 
black) of the AFM interface layer after the initial cooling procedure for dilution p = 0.3. AFM 
thickness t = 1 (upper figure) and t = 10 (lower figure). Vacancies are left white. 

FIG. 3. Exchange bias field versus thickness t of the AFM layer for different values of the 
dilution. 
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Figure 1: Nowak et al 
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Figure 2: Nowak et al 
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Figure 3: Nowak et al 
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